General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



Polytechnic 

institute 

d$ feylsrte 


Department of 
Aerospace Engineering 
and Applied Mechanics 

7 

/ 

/ 



J-wii-1*9o~/1j SfiiAfcAilttu *t> l TURfcUifcM N77-1*J30 

uOJAi.A*I L* it* CAi.CUi*IitRS A w i. I WC C A A I 

(fox/tkcni , ic lost, ot Ney 

loifc, i-a 4 . u. a 1,^04 it .) J 1 y bC AOj/flf AU 1 UhCiaj. 

CSC I *0l JJ/J* 2 1 o Jt> 


SEPARATING AND TURBULENT BOUNDARY LAYER 
CALCULATIONS USING POLYNOMIAL INTERPOLATION 


By S. G. Rubin and S. Rivera 


FEBRUARY 1977 




Grant No. NSG 1244 

Approved for public release; 
distribution unlimited. 

POLY-M/AE Report No. 77-4 


SEPARATING AND TURBULENT BOUNDARY LAYER 


CALCULATIONS USING POLYNOMIAL INTERPOLATION 

by 

S. G. Rubin and S, Rivera 


This research was sponsored by the National Aeronautics 
and Space Administration, Langley Research Center 
Hampton, Va., under Grant NSG-1244 


POLYTECHNIC INSTITUTE OF NEW YORK 
AERODYNAMICS LABORATORIES 
February 1977 


POLY-M/AE Report No. 77-4 


SEPARATING AND TURBULENT BOUNDARY LAYER 


CALCULATIONS USING POLYNOMIAL INTERPOLATION t 

by 

S. G. Rubin* and S. Rivera 

Polytechnic Institute of New York 
Aerodynamics Laboratories 
Farmingdale, New York 

ABSTRACT 

Higher-order numerical methods derived from polynomial 
spline interpolation or Hermitian differencing are applied 
to a separating laminar boundary layer, i.e., the Howarth 
problem, and the turbulent flat plate boundary layer flow. 
Preliminary results are presented herein; as in earlier non- 
separating laminar flow studies, it is found that accuracy 
equal to that of conventional second-order accurate finite- 
difference methods is achieved with many fewer mesh points 

a 

and with reduced computer storage and time requirements. 
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Introduction 

(1-4} 

In a series of recent reports ' ' polynomial spline and Bermite 

interpolation procedures have been developed and applied to laminar 
boundary layer and incompressible Navier-Stokes problems. The poly- 
nomial interpolation development is quite general so that a variety 
of second, fourth and sixth-order spatially accurate representations 
for both convection and diffusion terms are possible? e.g., a quadra- , 

tic polynomial leads to the familiar central finite-difference dis- 

( 2 ) 

cretization formulas . A detailed discussion of the various form- 
ulations with comparisons, interrelationships and examples is given 
in Refs. (1-4) , All of the governing systems are block tridiagonal; 

the unknowns at each mesh point <n. are at most the function V. and 

3 3 

the spline approximations m^, of the derivatives (V^)j, (V^j, 

respectively. For the second-order procedures, the governing system 
can be reduced to include only V\; for the fourth-order procedures, 
a reduction requiring only Vj and Mj is possible, and for the sixth- 
order development the complete 3x3 block tridiagonal system is 
required . 

In the present note, the applicability of the polynomial spline 
methods for a laminar separating flow, i.e., the linearly deceler- 
ating Howarth problem, and the turbulent boundary layer over a flat 
plate is investigated. In the latter case, the Cebeci-Smith two- 
layer eddy viscosity model is used for turbulent closure. For both 
the Howarth and turbulent flat plate flows the solutions have been 
studied extensively by many investigators so that meaningful compari- 
sons with the present calculations are possible. These two problems 

i 
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provide somewhat different and in certain ways more severe tests of 
the polynomial interpolation methods than did previous boundary 
layer calculations, see Refs. (2-4) . From these earlier laminar 
flow studies, it was found that the fourth-order spline methods 
generally required one-quarter the number of grid points, in each 
coordinate direction, as did second-order central finite-difference 
methods in order to achieve equal accuracy This can lead to 

reductions of from 50% to 75% in computer time and storage 
The present calculations will determine the gains, if any, when 
points of zero shear are approached and when a thick turbulent 
boundary layer with a thin high shear viscous sublayer is to be ap- 
proximated . 

Four of the polynomial interpolation procedures will be 

considered here: (1) Second-order central differences; (2) a 

mixed second and fourth-order method termed Spline 2, (3) a fourth- 
order system termed Spline 4, and (4) a sixth-order Hermite 6 de- 
velopment. The block tridiagonal governing equations are presented 
in the following sections. For details of the derivations, see 
Refs. (2-4) . 

Governing Equation s 

After appropriate normalization and transformation, the govern- 
ing equations for both the Howarth flow and turbulent flat plate 

( 5 ) 

boundary layer become v 


(1+e) V + ( cv+f+2 ?f .) V +8 ( 1-V 2 ) =2 |W e , 

ryn s <n s 


£ n- v ' 


(1) 

( 2 ) 


where V-V( n)=u/u e ; u is the streamwise velocity and e denotes po 


2 


tential flow conditions at the boundary layer edge 


£ = £(?, T,) , U e - U G (?) » «=* 


2?(uJ 


e' x 


u 


X 


S = 5(x) = J u £ (x)dx, r> = n(x,y) = 


y u . 




o - (2 v5) 

x, y denote the axial and normal coordinate directions; v is the 
kinematic viscosity; e is the eddy viscosity for the turbulent 


boundary layer. The two-layer Cebeci-Smith 


( 6 ) 


model is used for 


e, see Ref. 7. For the Howartli flow: e=o, u = 1-x = (1-2 f) % 

For Turbulent flow: u ^constant, 0=o. 

e 


P ~2 ?/ ( 1-2 


Polynomial Spline Interpolation 

Consider a mesh with nodal points at T>=r>j . Define the mesh 

width h. = ru-n. , and cr.=h . . .,/h . . At the mesh points, 

3 3 1-1 3 3 

V( 5/ n) == V( 5, rtj)=Vj . Define an n order polynomial S(?,'tl) on 
Tj-1, jl such that S ( ?, r»j) =Vj, and itk, denote the polynomial ap- 

proximations of the functional derivatives (V^) ., (V^) ^ , respective- 
ly. A similar procedure is applied on [j , j+1 3 and continuity of 
S(Sf'n) and its derivatives is prescribed at n j . This leads to a 
coupled system of equations for V^, nu, at each mesh point. The 
system is closed with the differential equation (1) for ; func- 
tional derivatives are replaced by their polynomial approximations . 

A similar procedure is used for the function f(§, n), Finally, in 
the marching or | direction, simple backward differencing is used 
for V ? ;i.e., (V j n+1 -V : . n ) /A|= (V^) ^ , where ?=nA?; n=0,l,2,. 

For the polynomials considered here, the governing block tri- 
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I 
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diagonal systems for m. , M., with o.t=iif=l, are as follows. The cleri* 

J *] J 

vations and results for tf^l are found in Refs. (1-4) . 

A j I, j+l +B jLj + CjLj.x* , D ; ., where L^rm.M]* , (3« 

Central Differencing: A^= [jj q],B_.= [ o 1 ] ' c j = [o 0.1 ( 31 

spline 2,4*= 6 Rj = [J °], 3B.= [* °], 6 Cj = [J °j (3, 

Hermite 6= 120a.- [f/ h _]£], 1!' [° *], 120C = [‘« /h $ 


. 2 n r (v j+r v j-i )h 

2h D.= J J 


3 2 ( V j+l- 2V 3 +V j-l ) 


in all cases. 


Equation (1) becomes 


/ •= . \ n+ 1„ n+ 1 , / i n+ i n*t 1 . -nt 1 , .. 2 . n*t I n » nt 1 nt 1 , . n= t a\ 

(1+e . ) M. +(c +f+2 ?£ J ra. + fl (1-V. ) =25.V. (V. - <?■ ) (4) 

jjnSjj 3 3 j j J 

aT 

2 1 * 

Non-linear terms (V^) 11 ' are treated by quasi-linearization; 

•! « 7 2. n+1 n„ n+1 n 2. n 

i.e./ (V.. ) ~2Vj Vj -(Vj ) . 

The block-tridiagonal system (3a) and (4) are solved by a now 
standard inversion algorithm ^ . For Spline 2 or 4 the blochs can 
be reduced to 2x2, while for Hermite 6 they are 3x3. The implicit 
formulation (4) is unconditionally stable. 


*For Spline 4, (V^) j is approximated by ^-2Mj)/l2, see 

Refs . (2,3). 
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j RO "Ults . 

^ The results for both the Howarth and turbulent flow examples 

are quite encouraging. Significantly fewer grid points are re- 

• quired to achieve accuracy comparable with that of finite-dif fer- 

, ence solutions. Even though operational counts are increased due 

to the larger matrix systems, reductions in computer storage and 

time are only , slightly less than that found in earlier laminar 

studies . 

Howarth Problem 

; Typical results are shown on Table 1 and Fig, 1. N denotes 
’ the number of mesh points and h ^= is the mesh width at the sur- 
face. Even with the very coarse grid (hj-1.0), the Spline 4 re- 
; suit is within 4% of the usually accepted value of x==0.1198^? 

- moreover, the Hermite 6 solution is in error by less than 1%. Al- 
so shown on Table 1 are the results ^ of an alternate but somewhat- 
less accurate fourth-order scheme termed here Hermite 4 or compact 
■ This scheme is discussed in Refs. (3, 4, 8) . The finite-dif- 
ference results are in error by more than 20% and require at least 
four times as many grid points to reduce this error to 4%„ It 
< would appear then that separation can be predicted accurately with 
a minimal number of grid points if polynomial interpolation methods 
are applied. 


N 

h 2 

a j 

F.D. 

.HERMITE 4 . 

SPLINE 4 

HERMITE 6 

' 61 

0.1 

1.0 

0.1199 

0.1198 

0.1198 

0.1198 

7 

1 

1.0 

1.0 

0.1458 

0 .1121 

0.1159* 

0.1193 


*If second-order boundary conditions are applied, this value is 0.1191. 
Since the shear vanishes near separation, and Vrmrvn ( §, o) ~Vri( P, o) , these 
boundary conditions are in fact fourth-order at separation. The spline 4 
boundary conditions can be modified .toQU.se V Tvnnrti; f, o) throughout and im- 
prove the accuracy of the solution^ ' This is important for coarse 

grids and near separation points. 

TABLE I SEPARATION POINT - HOWARTH PROBLEM 
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Turbulent Flow 


Similar trends are observed for the -.urbulent boundary layer, 
although considerably larger number of me3h points are required 
to approximate these thicker high .tfiear regions. Coordinate trans- 
formations might facilitate further reductions in the number of 
points. These were not considered here. Typical results are shown 
on Table 2 and Figs. 2,3. All indications are that reductions of 
from one-third (Spline 4) to one-quarter (Hermite 6) the number of 
mesh points, and 50% in computer time and storage, are possible. 


N 

h 2 

CT. 

3 

Cent. Diff* 

Spline, 4 Hermite 

141 

0.001 1.063 

3.456 

3.455 

61 

0.2 

1.0 

5.479 

3.891 3.476 

61 

0.1 

1.0 

5.148 

3.682 

21 

0.02 

l+4h . 
3 

4.396 

3.472 

11 

0.05 

1.5 

6.823 

3.759 

11 

0.05 

l+4h . 
3 

- 

3.598 

Wieyhardt 

Data 

(Ref. 9) 

C_ = 3.45 x 10 
f 

-3 




^ C f^ laminar 

0.664 x 10~ 3 

T&BLE 2: 

COEFFICIENT OF 

FRICTION (C £ x 

10 3 ) at R x = 10 6 , 
© 


At R g x = 10^ the turbulent flat plate boundary layer is approximate- 
ly twice as thick as its laminar counterpart and the local shear 
stress is five times greater. The Spline 4 results with 11 and 21 
points are quite good. Hermite 6 solutions are excellent but are 
given here only for *= 1,0; the variable grid formulations are given 
in Ref 0 3 and will be applied in future studies. In view of the 


6 



present: results, the application of higher-order polynomial interpo- v 
lation for separating and turbulent flows appears to be quite promis- 
ing. More detailed results and discussion can be found in Ref. 7. 

In the Appendix attached herein are a number of additional figures from 
Ref. 7 describing the flow field and solutions for both laminar and 
turbulent cases. Further analysis of turbulent flows with uniform 
and variable grid spline 4 and Hermite 6 methods are presented in 
Ref. 10. Several transformed coordinate systems and flows with surface 
mass transfer are also considered. 
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fig. 3 STREAM WISE VELOCITY VARIATION- 
TURBULENT BOUNDARY LAYER 
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